function [output] = fn_alphaN(R1,lambda,phi,s)

rho1_s = fn_rho1(s,phi);
rho2_s = fn_rho2(s);

output = rho2_s*(rho1_s - lambda*R1)/((1-lambda)*R1);

end